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We present estimates for the 3D Ising model on the cubic lattice, both regarding interface and bulk properties. 
We have results for the interface tension, in particular the amplitude ao in the critical law a = aot'^, and for 
the universal combination R- — a^^. Concerning the bulk properties, we estimate the specific heat universal 
amplitude ratio A+/A-, together with the exponent a, the nonsingular background of energy and specific heat 
at criticality, together with the exponent v. There are also results for the universal combination /s^^, where /a 
is the singular part of the free energy. Details can be found in [Q] (interface) and Q (bulk). 



1. INTRODUCTION 

In this conference contribution, we summarize 
recent efforts to determine various critical 
quantities of the 3D Ising universality class using 
the Monte Carlo (MC) method. These efforts are 
worthwhile for at least three reasons: 1) There 
are many systems in nature belonging to the Ising 
universality class. Thus a comparison of theoret- 
ical predictions with experimental results is pos- 
sible. 2) Getting precise estimates is still a chal- 
lenge, both from the algorithmic side and from a 
theoretical point of view. MC data with increas- 
ing statistical precision put standard assumptions 
to test and may require more and more refined 
theoretical modelling. 3) MC precision has now 
become competitive with, e.g., high/low temper- 
ature series expansions, the e-expansion and 3D 
field theoretic calculations. Comparing the re- 
sults is interesting and might point out sources of 
systematic errors in either method. 

2. INTERFACE PROPERTIES 

Properties of interfaces separating extended 
domains of different phases are of great interest, 



both from the experimental and the theoretical 
point of view. For recent experimental work mea- 
suring observables that are related to and can be 
compared with quantities studied in the present 
work, see 

Figure 1 

Lattice geometry 
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We consider the 3D Ising model on a cubic lat- 
tice of size LxLxT with spins Sx = ±1, Hamilto- 
nian H{s) = -(3 Yl,<x,y> ^^^yy where <xy> de- 
notes a pair of nearest neighbour sites on the lat- 
tice, and partition function Z — J2 s i^)] ■ 
By Zp we denote the partition function with peri- 
odic boundary conditions in all three directions. 
Antiperiodic boundary conditions in i-direction 
favour the occurence of interfaces (perpendicular 
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to the t-axis), see figure 1. The corresponding 
partition function is denoted by Zq. If we can 
assume that there is exactly one interface in the 
system, its free energy is given by 

F,^Fa-Fp + lnT, (1) 

where hiT takes care of the free motion of the 
interface in the i-direction. 

We obtained the interface free energies by inte- 
gration of the interface energy Eg = {H)a — {H)p: 

rP 

Fs = FsiPo)+ / dp'EsiP'). (2) 

Various state of the art algorithms were em- 
ployed: the boundary flip cluster algorithm for 
the determination of Fs{(3q), the WoliT single clus- 
ter, the Ising interface cluster and the multispin 
coded demon algorithm for the determination of 
Eg . Detailed descriptions and the appropriate ref- 
erences can be found in Q . 

Interface tensions a were then obtained by fit- 
ting the free energies with either 

Fs~Cs + cjL^ (3) 

or I 

FgCiCs + aL^ +l/{4aL'^). (4) 
We chose lattice sizes 323, 433^ 543 g^^^ gg3 jy^^^g^ 
were taken on a dense grid of f3 values in the range 
0.223 . . . 0.23. In table |] we quote a few of our fit 
results for the interface tension. Here, "fit2" la- 



Table 1 



Selection of a 


estimates 




fit 


/3 ^ 0.2240 


f3 = 0.2255 


fit2 


0.004784(9) 


0.008810(9) 


96vs64 


0.004799(17) 


0.008822(18) 


fit3 


0.004782(6) 


0.008801(6) 



bels a fit with eq. (||), discarding the L = 32 data, 
under the name of "96vs64" goes a determination 
of Cs and a from the L = 64 and L = 96 data 
alone, and "fit3" refers to a fit of data from all 
four lattice sizes with eq. (||). The results from 
the different fit types are fairly consistent. 



Our next task was to fit the interface tension 
results to the critical law 

a{P) -cTof^il + agt" + ait). (5) 

Here we included the leading corrections to scal- 
ing terms. We compared the results stemming 
from the two definitions of the reduced tempera- 
ture, ti = 13/ — 1 and t2 = 1 — P. Further- 
more, we fixed the following parameters: (3c = 
0.2216544(6) §], ^i = 1.262, and 9 = 0.51 §. 
These were the most precise estimates we could 
find in the literature. The fit results are quoted 
in table ||, again for three different fit schemes. 



Table 2 



Fit results for do 7 o-e and ai 



fit 


t-def 


cro 


ae 


ai 


fit2 


h 


1.549(11) 


-0.409(72) 


0.01(20) 






1.549(11) 


-0.397(76) 


1.13(21) 


fit3 


h 


1.5428(73) 


-0.376(49) 


-0.06(14) 






1.5421(73) 


-0.362(51) 


1.06(14) 


96vs64 


h 


1.571(20) 


-0.57(13) 


0.47(36) 






1.571(21) 


-0.56(13) 


1.61(37) 


Note 


that 


the (To 1 


results for 


the two 


t— definitions 


are always 


consistent 


with each 



other, ai should have a jump oi ^ k, 1.26 which 
is obeyed within the errors. Taking into account 
the various sources of uncertainty (including the 
variation with the type of fit) we arrive at the 
final estimate uo = 1.55(5). 

We finally employed our cr-data in an estima- 
tion of the universal ratio Including correc- 
tions to scaling, it can be defined through 

i?_ -a(/3)C(/3)2-cr", (6) 

where (3 ^ (3c from above, lo is the correction 
to scaling exponent. We fixed it to 0.81(5) 
The ^-data were taken from |6|. Our result is 
R- = 0.1040(8), where the error is mainly due to 
the uncertainty in the exponent uj. A compari- 
son with some results from the literature is given 
ini. 
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3. BULK PROPERTIES 

We consider a 3D Ising model on a cubic lattice 
of size L^, with periodic boundary conditions in 
all three directions. Important quantities are the 
free energy per volume 



-1/(3^3) In Z, 



(7) 



the energy, E — —df/d/3 and the specific heat, 
C = dE/dp. The quantities are split in a non- 
singular and a singular part, e.g., E — E^s + Eg. 
By standard arguments the following finite size 
scaling laws are expected to hold: 

E ~ £:„s + constE L-'^+i/'^ 
C ~ Cns + constc L-'^+^l" 



at criticality, and 

E ~ i?ns — Cns Pet T A± Pc 



I- a 



(8) 



(9) 



above and below the transition. Using state of the 
art algorithms we obtained estimates for energy 
and specific heat at (3c = 0.2216544 for L ranging 
from 12 to 112. The relative error of the energy 
is of order 10~^ throughout. We then did a com- 
bined fit of energy and specific heat to eq. (||). 
Our results for Ens and Cns are given in ta- 
ble ^ Lmin is the smallest lattice size used for 



Table 3 

Fit results from critical energy and specific heat 





X 


V 


Ens 


Cns 


12 


5.6 


0.6316(4) 


0.330229(5) 


-12.13(33) 


16 


2.25 


0.6315(8) 


0.330218(7) 


-11.83(61) 


20 


0.75 


0.6308(10) 


0.330209(8) 


-11.12(76) 



the fit, and X denotes per degree of freedom. 
We find the precision of the z/-estimate of the last 
line of the table quite remarkable. 

In the off critical case, we always aimed at the 
infinite volume limit of the energy. Our data 
cover the range from (i = 0.21971 to /3 = 0.2230. 
Lattice sizes up to L = 128 were employed. 
The results of fits with eq. (^ are given in ta- 
ble | A gain we used two definitions of i, namely 



Table 4 

Fit results from off critical energy 





h 


t2 


a 


0.1115(37) 


0.1047(48) 


A+IA^ 


0.550(12) 


0.567(16) 


Ens 


0.33029(3) 


0.33030(4) 


Cns 


-9.86(80) 


-11.37(1.17) 



ti = 1 — P/Pc and t2 = Pc/P—l- There is a slight 
mismatch of the estimate for ii^ns with the esti- 
mates from the critical data in table ^. It is most 
likely due to corrections to scaling that have not 
been taken into account. 

We finally estimated the universal constant 
fs on both sides of the transition. Estimates 
for fs were obtained by integration over P of the 
singular part of the energy, obtained from the ap- 
proximation Eg K E — Ens — Cns(/3 — Pc)- Us- 
ing our background estimates £^ns and Cns to- 
gether with the correlation lengths given in 
we obtained /s^^ = 0.0355(15) for P / P^ and 
/s^^ = 0.0085(2) for P\Pc. 

4. CONCLUSIONS 

For various interface and bulk quantities we 
could significantly improve on older numerical es- 
timates. Especially interesting is the quite precise 
determination of v from the finite size scaling of 
the energy and specific heat at criticality. 
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